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O : Abstract 

r-| ■ The Balitskii-Fadin-Kuraev-Lipatov (BFKL) and the Gribov-Lipatov-Dokshitzer-Altarelh-Parisi 

ji. . (GLDAP) evolution equations for the diffractive deep inelastic scattering at - ^ 1 are shown 

^ , to have a common solution in the weak coupling limit: F2{x,Q^) oc [asiQ"^)]^^ (-) • The 



exponent 7 and the pomeron intercept Ajp are related by 7A]p = | for the Nf = 3 active 



H ' flavors. The existence of this solution implies that there is no real clash between the BFKL and 

GLDAP descriptions at very small x. We present derivation of this solution in the framework 
of our generalized BFKL equation for the dipole cross section, discuss conditions for the onset 
of the universal scaling violations and analyse the pattern of transition from the conventional 
Double-Leading-Logarithm approximation for the GLDAP evolution to the BFKL evolution at 
large ^ 
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1 Introduction. 

The asymptotic behavior of diffractive scattering in perturbative QCD is usually discussed in 
the framework of the BFKL pomeron [1]. Recently, we took advantage of exact diagonalization 
of the S-matrix of diffractive scattering in terms of the dipole cross section [2,3] and derived the 
generalized BFKL equation for the perturbative dipole cross section [4-6]. In [4-6] we studied 
the spectrum of our generalized BFKL equation and determined the intercept of the pomeron 
and the asymptotic dipole cross section for the rightmost singularity in the complex j-plane. 

The central result of the present paper is a derivation of the asymptotic pomeron solution of 
our generalized BFKL equation. 



(^Tp{^,r) = crip(r)exp(Aip^) oc r^ 



7-1 



exp(AipO , (1) 



_as[r) 
where the exponent 7 is related to the pomeron intercept Aip by 

Here /3o = H — ^Nf = 9 for Nf = 3 active flavours, r is the transverse size of the colour dipole, 
as{r) is the running QCD coupling and C, = log(— ), where Xq ~0. 1-0.01 corresponds to the 
onset of the leading-log(-) approximation. We shall demonstrate that, to an accuracy ~ ^^(r), 
the new solution (|l]) is a low-x limit of solutions of the conventional GLDAP equation [7] . This 
new solution must be contrasted to: 

i) The Double-Leading-Logarithm Approximation (DLLA) solution [4] to the GLDAP evo- 
lution equation. 



n ■ - 

(^dlla{^, r) = a(0, r) ^ — — — - ~ r^as{r) log 



where 77 is the expansion parameter of DLLA, t] = ^^L{r) and L{r) = log 
ii) The solution of the scaling BFKL equation [1] in the case of fixed as, 



^exp(2v/^), (3) 



as(r) 



crjpitr) oc rexp(AipO • (4) 

The conventional DLLA solution (^ sunims the leading powers ^"L(r)", our new solution (|1]) 
manifestly summs all powers of L{r). 



The long going debate in the hterature on the clash between, and transition from, the GLDAP 
to the BFKL regime with rising -, is centered on a comparison of the DLL A solution (|^) and the 
scaling BFKL solution (^ (for the recent review and references see [8,9]). The existence of our 
new solution of the generalized BFKL equation has the two major implications: Firstly, there 
is no real clash between the GLDAP and BFKL evolutions and the GLDAP evolution remains 
a viable description of deep inelastic scatering at very large ^ in the perturbative regime of 
c^siQ"^) *^ 1- Secondly, comparison of solutions (|l]) and (§) shows the dramatic difference 
between the cases of the running and fixed strong coupling for the small-x behavior (see also 
below, a discussion of the so-called DLLA identity). We conclude this introduction citing the 
universal structure function 

1 V /w^^ 

F,{x,Q')^ — — - , (5) 

which follows from the dipole cross section (|l|). Most unfortunately, the onset of this universal 
regime only is expected beyond the kinematical range of the HERA experiments. 

The paper is organized as follows: We begin with the brief review of the lightcone description 
of deep inelastic scattering in terms of the dipole cross section and present our generalized BFKL 
equation for the dipole cross section. Then we discuss a transition from from the DLLA solution 
to GLDAP equation to the BFKL solution. We introduce a new comparison of the GLDAP 
and BFKL solutions in terms of the so-called DLLA identity, and present a derivation of our 
new solution (|1]) and of the universal scaling violations at asymptotically large -. Finally, we 
comment on the significance of the consistent use of the running coupling for the pomeron cross 
section. In the conclusions section we summarize our main results. 

In the region of very large - the virtual photoabsorption can be viewed as interaction with 
the target proton (nucleus) of the multipartonic lightcone Fock states (gg, qqg.--) of the photon, 
which are formed at large distance Az ~ -^— upstream the target. Consequently, the longi- 
tudinal momentum partitions Zi and the transverse separations pi of partons are conserved in 
the scattering process, and the diffractive S-matrix is exactly diagonalized in the (p, z) repre- 
sentation [2]. Photons do not couple to gluons and the higher qqgi...gn states are radiatively 
generated from the qq dipole of size r. In [4] we gave regular procedure for construction of the 
corresponding multiparton lightcone wave functions and of the multiparton total cross sections. 



The Fock states with n gluons give the oc log'^(-) contribution to the total photoabsorption cross 
section, which can be reabsorbed into the energy dependent dipole cross section 

^i^,r) = Y.^cT^ir)C. (6) 

The total photoabsorption cross section can be written as an expectation value [2] 

aTAl*N,^,Q') = J^ dz J d^r\^TA^,r)\M^,r) (7) 

over the wave function of the qq Fock state. The wave functions of the (T) transverse and (L) 
longitudinal virtual photon of virtuality Q^ were derived in [2] and read 

\^T{z,r)\' = ^E^Kt-' + (1 - ^WK,{err + mjMerr} , (8) 

\^L{z,r)\^ = ^E4^? Q'^'i^ - zfK,{erf , (9) 

where Ki{x) is the modified Bessel function, e^ = z{l — z)Q'^ + m'j, ruf is the quark mass and 
z is the fraction of photon's light-cone momentum carried by one of the quarks of the qq pair 
{0 < z < 1). Then, the structure function is calculated as F2(^,Q^) = Q^ [ax + o-^] /{An^aem)- 
Notice that the dipole cross section <j{C,,r) is universal, only I^E't.lP depend on Q'^ and m'j. 
In [4-6] we have shown that o"(^,r) satisfies the generalized BFKL equation 

^^|l^ = /C®a(e,r), (10) 

where in terms of the expansion (|^) the kernel /C is defined by 



'^ f j2^ 2 
— Jdp.fia 



(T„+i(r) = ]C®(Tn{r 

2 

gsiRi)Ki{fiGpi)— - gs{R2)Ki{pGp2) — 

pi P2 



0-„(pi) + 0-„(p2) - (yn{r)] . (11) 



Here Re = l/pc is the correlation radius for perturbative gluons, Ri = min{r, pj}, gs{r) is the 
effective colour charge, 

where C ~ 1.5 [2] and we impose the infrared freezing as{r > Rj) = ag = 0.8 [6]. 



In the BFKL scaling limit of r, pi ^ Re and of fixed as, the kernel /C becomes indepen- 
dent of Re, Eq. (|T0|) becomes equivalent to the BFKL equation [1] and has the familiar BFKL 
eigenf unctions 

E{^,r) = r'+''^exp[^A{tu)] (13) 

with the BFKL eigenvalue (intercept) 

AM = ^[2M/(i) - nl - ^^) - nl + ^^)] , (14) 

TV II 

where '^[x) is the digamma function. 

Deep inelastic scattering at large Q^ probes a[^, r) at small r^ oc 1/Q^ ([2-4], see also below). 
In the DLLA of large but finite ^ and large lj{r\ the kernel /C is dominated by r^ <^ pj <^ R^f, 
one can neglect a{r) ^ cr(pj) and factor out as{r) in Eq. (|lT]), which in the DLLA takes on a 
particularly simple form [4-6] 

a„+i(r) = /C ® cxn{r) = / -^an{pj , (15) 

TC Jr^ p 

what is equivalent to the large-- limit of the GLDAP evolution equation [7]. As a boundary 
conditon, one can start with the dipole cross section for interaction with the nucleon target [2] 

Mr) = y y j^^^^j—as{e)as{^')[l - G,{3k')] [l - exp(-^fcf)] , (16) 

where Gp{q^) is the charge form factor of the proton, h? = iaax{k'^,^}. Then, iterations of 
Eq. (|1^) produce the familar DLLA solution (^ [4]. 

The question of whether there is a strong, experimentally observable, difference between 
the BFKL and GLDAP evolutions (the latter usually being considered to DLLA), is being 
discussed in the literature for quite a time ([8,9] and references therein). A comparison of 
the DLLA iterations (Jn+i{r) = /C o"„(r) = 12L(r)cr„(r)/(ri + l)f3o with the BFKL iterations 
<^n+i{r) = K ® anir) = Aip(T„(r) suggests the DLLA breaking at L{r)/n ^ L{r)/ y/fj < |Aip. 
(At large rj the dominant contribution to the DLLA cross section (^) comes from n ~ ^/Tj.) We 
can go one step further and compare in detail the ^-dependence of the solution of our BFKL 
equation and of the DLLA solution of the GLDAP equation, starting with the identical initial 
condition Eq. (|l^). (For the sake of definitness we consider Re = 0.275f which gives Aip = 0.4 
[6].) Specifically, we compare the effective intercepts Ae//(^,r) = d\oga{C,,r)/dC, for the two 



solutions. Since the DLL A asymptotics (|^) works at ,^ > 1, we shall consider the cross section 
( [T^ ) as a result of the GLDAP evolution by ^o units from a lower energy, i.e., we take for the 
DLLA solution 



crDLLA{^,r) = (Toir) 



eo 



■exp 



2JlHr) 



eo + e 



(17) 



eo + e 
It satisfies cr£,LLA{C, = 0, r) = (jQ[r) by the construction, and gives the DLLA effective intercept 

1 



A 



DDLA 



(e,r) 



4L r) 



(18) 



^3(^0 + 2(^0 + 0' 
The same cro(r) is taken as the boundary condition for the BFKL equation (jT^). (Here we 
assume it to correspond to x = Xq ~ 3 ■ 10~^, the more detailed BFKL phenomenology of deep 
inelastic scattering will be presented elsewhere [10]). 

Firstly we compute Ae//(^ = 0, r) from our BFKL equation (|1^) and make the readjustment 



Lir) 



log 



a 



ifr) 
S 



as[r) 



+ c 



(19) 



such that AciLLAi^ = 0, r) of Eq. (|TB|) gives a good approximation to this Ae//(^ = 0, r) at small 
r. (Recall that L{r) is defined up to an additive constant c < 1.) With ^o = 1-25 this is achieved 
taking c ~ 0.05. Secondly, we study how the BFKL and DLLA effective intercepts diverge at 
large ^ (Fig.l). At ^ = log(— ) ~ 1, both the DLLA and BFKL effective intercepts are smaller 
than Aip at r > 0.2f and are larger than Aip = 0.4 at smaller r. The good matching of the 
BFKL and DLLA effective intercepts at small r is not surprizing, since our generalized BFKL 
equation (|T0| , |TT1) has the GLDAP equation as a limiting case at small r [4-6]. With rising -, the 
BFKL effective intercept flattens and tends to Ajp, rising at large r and decreasing at small r, 
whereas the DLLA intercept monotonically decreases with ^ at all r, until the DLLA breaking 
Adlla(^,^) < Aip takes place at ^ = ^c{r) given by 



^c{r) = log 






3A2, 



log 



<ys{r) 



(20) 



(The much discussed boundary suggested in [11] is not born out by our accurate comparison of 
the DLLA and BFKL solutions.) The intercept Aip is small, Aip = 0.4 at fic = 0.75GeV [6]. 
The resulting large numerical factor 



which emerges in the r.h.s. of Eq. pO|), explains why the close similarity of the BFKL and DLLA 
solutions persists in such a broad range of r and x, relevant to the HERA experiments. 

The dipole cross section cr(^, r) is related to the more familiar gluon structure function G(^, r) 
by [4,12] 

Gi^,r)=xgix,r) = ^^^ (22) 

where g{x,r) is the density of gluons at x = xoexp(— ^) and the virtuality Q^ ~ 1/r^. In terms 



of the G(^, r) the DLLA equation (|T5|) is equivalent to 



which we shall refer to as the DLLA identity. One can easily evaluate k,{^, r) for the experimen- 
tally measured gluon distributions. It is interesting to look at the possible departure from the 
DLLA identity of the above described solution of our generalized BFKL equation ([TO| , p]TD subject 
to the boundary condition (|T6|) . The results of such a test are shown in Fig.2. We find that our 
BFKL solution produces k{^, r) fti 1 in a very broad range of ^ and r of the practical interest: the 
DLLA identity holds to the few per cent accuracy at r < \Rc, and to the (20 — 30)% accuracy 
even at large r, up to r^ < |_R^. The somewhat oscillatory r-dependence of the Ae//(^,r) is 
quite natural and has its origin in the contribution of oscillating harmonics with large |i/| (for 
instance, see Eq. (0)), which die out at large C,. 

This remarkable finding of k,{C,, r) ^ 1 can be understood as follows: In the weak coupling 
limit we can factor out as{r) from the kernel /C, and the generalized BFKL equation (|1^,0) for 
G{C,,r) takes the form 



dG{^,r) 3 
= —^ I d Pi 



as{pi)G{^,Pi) , as{p2)G{^,p2) r^as{r)G{^,r) 



+ 



pI pI pIpI 



(24) 



di 27r2 

in which the leading contribution comes from p? ^ ^^- The major function of the term ex 
a5'(r)G'(^, r) is a regularization of the logarithmic singularity at pi^2 -^ 0, and at asir) ^ 1 this 
term can be neglected at the expense of the integration cutoff Pi2 ^ ^^- As a result, Eq. 



takes the form which is identical to the GLDAP equation (p!5|): 



dG{i,r) 12 rm 
oi Po Jo 



Now we wish to emphasize that making DLLA is not imperative when solving the GLDAP 
equation, and evidently Eq. ( pSj ) has the one-parametric family of the small-r eigenf unctions 

1 



^(7, ^, r) = exp[7L(r)] exp(AO 



17 



exp(AO , (26) 



_as{r)_ 

where the exponent 7 is related to the intercept A by Eq. (|2|). The corresponding dipole cross 
section is given by Eq. (|I|). In contrast to the conventional DLLA solution ( PB| ) which only 
summs the leading powers [L(r)^]", our new solution (^,0) manifestly summs all powers L{r)^. 
One can easily check that the neglected oc as{r)G{S,,r) term in Eq. ( P^ gives the oc as{r) 
correction to the solution (]26|). The small-r considerations alone do not fix the intercept A and 
the exponent 7, they are determined from the matching the solution (|ID with the large-r solution 
of our generalized BFKL equation. As we have shown in [6], the value of Aip is predominatly 
controlled by the semiperturbative region of r ~ i?c- 

In Fig. 3a we show the pomeron dipole section found in [6] by a numerical solution of our 
generalized BFKL equation. In Fig. 3b we show that these solutions to a good accuracy satisfy 
the property 

X = ^^ois{rV-^ = const. (27) 

Typically, this property holds up to r < \Rc- The smaller is the gluon correlation radius Re, the 
sooner starts, and the larger becomes with the increasing r, the oc 05 (r) correction to Eq. (|l|) 
(we keep a^ = 0.8). 

When solving the GLDAP equations exactly (numerically), one of course implicitly summs all 
powers of \og\i/ a s{Q'^y\. The obvious conclusion from the above derivation of (|26|,p]) is that the 
DLLA solution (^), which only is valid at moderate values of ^, evolves at larger ^, beyond the 
boundary (|2QD, into our new solution (|l]). Consequently, from the point of view of the practical 
phenomenology, there is no real clash between the GLDAP and BFKL evolutions. The BFKL 
evolution is evolution in ^ and requires as the boundary condition the knowledge of a{S, = 0,r) 
for all r at fixed ^ = 0. The GLDAP evolution is evolution in L{r), and requires the knowledge 
of a{C,, To) for all C, at fixed r = tq. Inspection of Fig. 1 shows that at 

r = ro ~ (0.3 - 0.7)i?, , (28) 

the DLLA and the BFKL intercepts are very close to each other and to the Aip in a broad 



range of x of the interest for the HERA experiments. Our conclusion is that choosing for the 
GLDAP evolution the boundary condition at tq ~ 0.15f, i.e., at Ql ~ 10-20GeV^ (for the Q'^ — r'^ 
relationship see below), 

G{x, Ql), xq{x, Ql), xq{x, Ql) oc x^ , (29) 

one will obtain the GLDAP solutions, which for all the practical purposes will be indistinguish- 
able from the BFKL solution, if A is taken equal to Ajp. 

Making use of properties of the modified Bessel functions, after the 2;-integration one can 
write 



ari^N, e, Q2) = 1^' dz J d'f \^t{z, r) | V(e, r) 



1 f^Hdr^ ai^,r) 1 rl°s"^(«'^w ^ ^., ^ .^n^ 

" Q-^ A/g. -^^^""q-^L ^"^^^^ ^^^' P^ ^'°) 

Notice, that the factor l/Q"^ in Eq. (^), which provides the Bjorken scaling, comes from the 

probability of having the qq fluctuation of the highly virtual photon. Making use of the pomeron 

soluiton (^) in (|30D, we easily obtain the universal asymptotic structure function (||). Since 

G{C,,r) slowly rises towards small r, the structure function F2{x,Q'^) receives a substantial (but 

not quite dominant) contribution from 

The explicit form of the wave function (H) leads to a rather large numerical factor i? ~ 10 in 

the relationship (^Tj) (for the related slow onset of the short-distance dominance also see [13]). 

According to Fig. 1, the BFKL effects are most significant at r > |i?c ~ 0.15f, and the scaling 

violations in the region of Q^ < (20-30) GeV^ seem to be the most interesting ones from the point 

of view of testing the onset of the BFKL regime. 

Eq. d^) predicts the universal, x-independent, scaling violation at small x: 

_ d\ogF2{x,Q^) ^ 4n dlog F2{x,Q^) ^ 

dlogasiQ^) PoasiQ') dlogQ^ ^^ ^ ^ 

The exponent 7 and the pomeron intercept Aip are not calculable within the perturbative 
QCD. For the guidance, we cite the results of an analysis [6]: 7 = 2.19, 2.81, 3.33, 3.70 for 
fiG = 0.3, 0.5, 0.75, 1.0 GeV, respectively. It would have been of great interest to see the uni- 
versal scaling violation (R2) experimentally, but closer inspection of Fig. 1 makes us to conclude 



that the apphcabihty region of ( p2D hes somewhat beyond the kinematical range of the HERA 
experiments. Indeed, the universal scahng violation comes along with the r-independent uni- 
versal x-dependence of the structure function (^), i.e, with Ae//(^,r) = Aip, whereas Fig. 1 
shows that in the kinematic range of HERA, the Aeff{C,,r) for our BFKL solution exhibits a 
still substantial r-dependence. 

The discussion of the BFKL effects in the current literature concentrates upon the approx- 
imation of fixed as (for the review and references see [8]). Our finding is that the effect of 
the running coupling constant is quite substantial. Firstly, the pomeron cross section (|l]) dra- 
matically differs from the oc r^ BFKL scaling solution (^,|13|) for the fixed as- Secondly, one 
can define the counterpart of the DLLA identity for fixed as too, making the straightforward 
substitution 



m-^_L ^«.(p)^£«.iog^ . (33) 



Then, making use of the BFKL formula for the pomeron intercept [1] 



Aip = A(0) = ^^^a^ , (34) 

one readily finds k = 2 log 2 for the BFKL pomeron solution (^. This departure from k = 1 
for the case of the running as{r) also emphasizes a dramatic difference between the cases of the 
fixed and running strong coupling. Therefore, the fixed-a;5 considerations are not appropriate 
for the phenomenology of deep-inelastic scattering. 

2 Conclusions 

The purpose of this paper has been a comparison of the BFKL and GLDAP evolutions at large - 
in the framework of the dipole cross section description of deep inelastic scattering. We found a 
new solution ([l|) which is common to the BFKL equation with the running QCD coupling and to 
the GLDAP equation, and provides a smooth matching of our generalized BFKL evolution and 
of the GLDAP evolution at small x. We derived the asymptotic form of the structure function 
(|^) and the universal law Eq. ( p2D for the scaling violations. Our principle conclusion is that 
the GLDAP evolution remains a viable phenomenology of the scaling violations at HERA and 
beyond, provided that one starts with the boundary condition ( pOj) at Ql ~(10-20) GeV^. 
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Figure captions: 

Fig.l - Comparison of effective intercepts of the DLLA solution (0) and of the solution of our 
BFKL equation. Both solutions start with the identical dipole cross section at x = 3 ■ 10"^. 
The pomeron intercept Ajp = 0.4 is shown by the horizontal line. 

Fig. 2 - Test of the DLLA identity for the solution of our generalized BFKL equation {Re = 0.27f, 
;xg = 0.75GeV). 

Fig. 3 - The left box: the pomeron dipole cross section ajp{r) for different values of fic- The 
straight lines show the r^ and r^ behavior. 
The right box: test of the scaling law x = ■\o"!p{i^)[(^s{f^)]'^^^ =const. 
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